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Measurements of event shapes and azimuthal correlations are presented for events 
where a Z boson is produced in association with jets in proton-proton collisions. The 
data collected with the CMS detector at the CERN LHC at ^/s = 7 TeV correspond to 
an integrated luminosity of 5.0 fb^^. The analysis provides a test of predictions from 
perturbative QCD for a process that represents a substantial background to many 
physics channels. Results are presented as a function of jet multiplicity, for inclusive 
Z boson production and for Z bosons with transverse momenta greater than 150 GeV, 
and compared to predictions from Monte Carlo event generators that include leading- 
order multiparton matrix-element (with up to four hard partons in the final state) and 
next-to-leading-order simulations of Z + 1-jet events. The experimental results are 
corrected for detector effects, and can be compared directly with other QCD models. 



Submitted to Physics Letters B 



*See Appendix A for the list of collaboration members 



1 



1 Introduction 

A detailed study of the production of a Z boson in association with jets in pp collisions at the 
CERN Large Hadron Collider (LHC) is of great interest. Measurements of this process can 
be confronted with the predictions of perturbative quantum chromodynamics (QCD) at the 
highest accessible energies and for a broad range of kinematic configurations. Considerable 
theoretical progress has been made in this field, such as developments in next-to-leading-order 
(NLO) calculations for up to four hard partons produced in association with a Z boson |1 J, NLO 
predictions for Z + 1-jet production that can be interfaced to parton shower (PS) approxima- 
tions |j2H6j|, and leading-order (LO) multiparton matrix-element (ME) event generators such as 
ALPGEN 0, MadGraph 0, and SHERPA \9\, with provision for PS development. In addition, 
Z + jets production corresponds to a major background to many other processes at the LHC, 
such as the production of top quarks, and it is important in searches for supersymmetric par- 
ticles and Higgs boson physics. An improved understanding of Z + jets production over the 
largest possible regions of phase space can therefore provide a helpful tool for extracting small 
signals. 

Previous studies of angular correlations between the Z and the "leading" jet (the one with the 
largest pj) and between the two jets of largest pj have been reported at the Tevatron by the 
DO Collaboration [10| and at the LHC by the ATLAS Collaboration using 36 pb^^ of integrated 
luminosity [llj. In this Letter, the comparison of models with data for highly boosted Z bosons 
with pj > 150 GeV is of particular interest. This region of phase space is critical in searches for 
new phenomena that are based on a large apparent imbalance in the total transverse momen- 
tum. Such imbalance can be produced, e.g., by the Z — > vt/ standard model (SM) background. 
The uncertainty of this background contribution is limited by the accuracy of current Monte 
Carlo (MC) models, which can be improved through studies of leptonic {£^£^) decays of Z 
bosons and their correlations with the associated jets. 

In addition to azimuthal distributions, we provide the first measurements of variables that 
categorize the topological structure of Z + jets events. Multijet production at e+e^ and ep col- 
liders was used in the past to tune parton showers and fragmentation functions in MC event 
generators, as well as to measure the values of the strong coupling constant l[T2l - [T6| . A set of 
event-shape variables suitable for hadron colliders has been proposed in Ref. I(T7[, which pro- 
vides resummed perturbative predictions at next-to-leading-log (NLL) for these variables. A 
measurement of event shapes in multijet events was reported recently by the Compact Muon 
Solenoid (CMS) Collaboration IITSll . 

This Letter extends measurements of angular correlations and event shapes in Z + jets events 
by probing the features of final states containing Z ^ £^i^ decays, where i = fiore. Such final 
states, often referred to as Drell-Yan (DY), include contributions from 7* and Z/ 7* interference 
terms arising from the irreducible background of virtual photons (7*) from qq — ?► 7* — )■ i^£^ 
processes. The data were collected with the CMS detector at a center-of-mass energy of 7 TeV, 
and correspond to an integrated luminosity of 5.0 fb ^ . The observed angular distributions and 
event shapes in Z + jets production are compared with predictions from several MC generators, 
and comprise the first study of this kind to be reported at the LHC. 

2 CMS detector 

The origin of the CMS coordinate system is chosen at the center of the detector, with the z axis 
pointing along the direction of the counterclockwise proton beam. The azimuthal angle is de- 
noted as (p, the polar angle as 6, and the pseudorapidity is defined as jj = — In [tan {9/2)]. The 
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central feature of the CMS detector is a superconducting solenoid of 6 m internal diameter that 
produces an axial magnetic field of 3.8 T. A silicon pixel and strip tracker, a lead tungstate crys- 
tal electromagnetic calorimeter (ECAL), and a brass/ plastic-scintillator hadronic calorimeter 
(HCAL) are positioned within the field volume. Iron and quartz-fiber hadronic calorimeters 
are located outside the magnetic field volume, within each endcap region of the CMS detector, 
at 3 < \t]\ < 5. Muons are measured using gas-ionization detectors embedded in the flux- 
return yoke outside of the solenoid. A detailed description of the CMS detector can be found 
in Ref. Ull. 

3 Monte Carlo simulation 

All production processes of concern, namely the Z + jets signal and backgrounds correspond- 
ing to top-antitop quark pairs (tt), dibosons (WZ, ZZ, WW), and W + jets events are generated 
with MadGraph (version 5.1.1.0), which provides up to four-parton final states and is in- 
terfaced to PYTHIA (version 6.4.24) [20] using the Z2 tune f2T] to implement showering and 
hadronization of the partons. The CTEQ6L1 [22] parton distribution functions (PDF) are cho- 
sen for these calculations. Alternative models for signal include (i) SHERPA (version 1.3.1) @ 
(with up to four-parton final states) using the CTEQ6m PDF l|22|l and the default tune, (ii) 
POWHEG [2-5 1 for generating Z + 1-jet events at NLO using the CTIO PDF [23 1 and interfaced 
to PYTHIA (version 6.4.24) with the Z2 tune for parton showering and hadronization, and (iii) 
stand-alone PYTHIA (version 6.4.24) with the Z2 tune. The cross section for the signal is normal- 
ized to match the next-to-next-to-leading-order (NNLO) prediction for inclusive Z production 
obtained with FEWZ fl^ and the CTEQ6m PDF f22]. The tt cross section is normalized to the 
next-to-next- to-leading-log (NNLL) calculation from Ref. [25]. 

The detector response is simulated using a detailed description of the CMS detector based on 
the Geant4 package [26], and the MC simulated events are reconstructed using the same pro- 
cedures used for data. During the data taking, an average of ten minimum-bias interactions oc- 
curred in each bunch crossing (pileup). The prevailing beam conditions are taken into account 
by reweighting the MC simulation to match the spectrum of pileup interactions observed in 
data. 

4 Event selection and reconstruction 

Event selection starts by requiring two high-px leptons at the trigger level. For muons, this cor- 
responds to an online pj threshold of 13GeV (17GeV during periods of higher instantaneous 
luminosity) for the muon of largest px (leading muon), and 8 GeV for the subleading muon. 
For electrons, the corresponding trigger thresholds are 17 GeV and 8 GeV. Offline, muon can- 
didates are reconstructed through a simultaneous fit to the hits recorded in the tracker and the 
muon detectors | |27| . Electrons are reconstructed using both calorimeter and tracking informa- 
tion Il28li . The two leptons of largest px (i-e., the two leading leptons) in the event are required 
to be of opposite electric charge and have px > 20 GeV, \rj\ < 2.4, and invariant mass satisfy- 
ing 71 < m££ < 111 GeV to be considered Z boson candidates. The lepton candidates are also 
required to be isolated from other energy depositions in the event. In particular, an isolation 
variable is computed using the scalar sum of transverse momenta of tracks and calorimetric 
energy depositions within a cone defined by AR = ^ (A^)^ + (A;^)^ = 0.3 around the di- 
rection of the lepton, where A^ is in radians. The contribution from pileup to this px sum is 
estimated from the distribution of the energy per unit area in the plane in each event [29], 
and is subtracted from the calculated sum. This corrected sum is required to be less than 15% 
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of the measured px of the lepton. Lepton reconstruction efficiencies are determined using sim- 
ulation, and corrected for differences between data and simulation using the "tag-and-probe" 
technique described in Ref. |30J . 

The inputs to the CMS jet clustering algorithm are the four-momentum vectors of the particles 
reconstructed using the particle-flow (PF) technique ||3TII32J , which combines information from 
different subdetectors. Jets are reconstructed using the anti-Zcj clustering algorithm ||33l , with a 
size parameter of R = 0.5, by summing the four-momenta of individual PF particles according 
to the FastJet package of Refs. Ii34ti35i . 

The reconstructed PF candidates are calibrated to account for any nonlinear or nonuniform 
response of the CMS calorimetric system to neutral hadrons. Charged hadrons and photons 
are sufficiently well-measured in the tracker and in the ECAL, and do not need such correc- 
tions. However, the resulting jets require small additional energy adjustments, mostly from 
thresholds set on reconstructed tracks and from the clustering procedure in the PF algorithm, 
but also from biases generated through inefficiencies in reconstruction. Jet-energy corrections 
are obtained using simulated events that are generated with PYTHIA (version 6.4.22), processed 
through a CMS detector simulation based on GEANT4, and then combined with measurements 
of exclusive two-jet and photon + jet events from data 1136 1 . By design, the jet-energy correc- 
tions correct reconstructed jets to the particle level II37I , as opposed to the parton level. An 
offset correction is also applied to account for the extra energy clustered in jets from the pres- 
ence of additional proton-proton interactions (in-time or out-of-time pileup) within the same or 
neighboring bunch crossings. The overall jet-energy corrections depend on the f] and pj values 
of jets, and are applied as multiplicative factors to the four-momentum vector of each jet. These 
factors range between 1.0 and 1.2, and are approximately uniform in rj. The jets accepted for 
analysis are required to satisfy pj > 50GeV and \ < 2.5. In addition, all jet axes are required 
to be separated by AR > 0.4 from those of lepton candidates from Z — )■ £^£^ decays. From MC 
studies, it is found that the selection efficiency of Z + jets candidates is almost independent of 
jet multiplicity. 



5 Observable quantities 

The observable quantities used to describe the properties of Z + jets events are the differential 
cross sections as functions of the azimuthal angles A(p{Z,ji) between the transverse-momentum 
vectors of the Z boson and the z* leading jet in the event; the azimuthal angles among the three 
jets of leading pj Acpdj,]]^), with / < k, and i and k corresponding to 1,2, or 3; and the transverse 
thrust Tj, defined as [17j 



Tt = 1 - max — ^ (1) 

where the four-momenta of the Z boson and the jets are used as inputs to calculate tj, with pT,i 
being the transverse-momentum vector of object /, and the sum running over the Z and each 
accepted jet in the event. The unit vector that maximizes the sum, and thereby minimizes 
Tj, is called the thrust axis. In the limit of the production of back-to-back Z + 1-jet events, 
Tj tends to zero (Fig. [l^). With additional jet emission (i.e., the appearance of a second jet), 
the values of thrust increase. Thrust is most sensitive to specifics of modeling of two-jet and 
three-jet topologies, while it is less sensitive to QCD modeling of larger jet multiplicities. For 
clarity of presentation, we display results in terms of In tj rather than tj. The largest possible 
value is reached in the limit of a spherical, isotropically distributed event, where In tj — ?■ ln(l — 
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II Ti) ~ —1 (Fig.[Tp). In this limiting case the term 2/ n originates from the uniform azimuthal 
distribution of the transverse momenta. 




Figure 1: Topology of Z + jets events: (a) for In tj — > — oo and A^(Z,;i) — > tt; (b) for In Tj — ?■ — 1 
and A(p{Z,ji) ^ n. 

To investigate the dependence of the topological properties on the complexity of the final state, 
the events are categorized as a function of jet multiplicity. In particular, the azimuthal distri- 
butions are reported in inclusive bins of one, two, or three jets. Furthermore, the phase space 
is characterized according to the pj of the Z boson, and measurements performed either for all 
Pj or in the region of pj > 150 GeV. Figure |2] shows distributions of the associated jet multi- 
plicity and the pj of systems for > 1 jet events, prior to background subtraction. Both 
distributions are presented at the detector level, and are within statistical uncertainty of the MC 
predictions for DY+jets, tt, and other electroweak (EW) background sources. It should be noted 
that pj refers to the transverse momentum of the Z boson, following background subtraction 
and unfolding of detector effects. 

6 Analysis procedure and systematic uncertainties 

The analysis procedure consists of the following steps: Z — > ee and Z ^ fiji candidates are 
selected as described in Section|4| the background is then subtracted and the resulting distribu- 
tions are unfolded to the particle level; finally, the two channels are combined. The dominant 
sources of systematic uncertainty arise from uncertainties in jet-energy scale, resolution of jet 
Pj, background subtraction, and the unfolding procedure. Individual steps of the analysis pro- 
cedure are detailed below. 

The Z + jets candidates include several sources of SM background (Fig.|2]), which are subtracted 
using predictions of the MadGraph MC event generator. The dominant background is from 
tt production, and is about 1.1%, 4%, and 8%, for the Njets > 1/ > 2, and > 3 inclusive bins of 
jet multiplicity. An independent evaluation of the background from tt events is also obtained 
from an eji control sample in data, which is selected by requiring the presence of an electron 
and a muon of opposite charge, but otherwise using the same criteria as used for selecting 
the Z + jets events. For each jet-multiplicity bin, the estimates obtained from data and MC 
simulation agree within 6%. The two estimates are consistent, given the uncertainties on the 
integrated luminosity (2.2%) L381 and on the tt cross section (6%) I25J . Other backgrounds 
(dibosons and W+jets) are much smaller, and are evaluated using MC simulation. A total 
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CMS, = 7 TeV L = 5.0 fb"' CMS, = 7 TeV, L = 5.0 fb"' 




Figure 2: Distributions for Z — )■ f/^w candidate events in data, compared with expectations from 
simulated signal and background contributions using MadGraph simulations normalized to 
the integrated luminosity of the data: (a) as a function of associated jet multiplicity Njets, and 
(b) as a function of pj of the dimuon pair (pj^') for Njets > 1- The dibosons WW, WZ, ZZ and 
W + jets backgrounds are collectively denoted as EW in the legends. The plots in (c) and (d) 
show the ratios of the data to predictions from MC. The error bars on the data points represent 
their statistical uncertainties. 

uncertainty of 10% is assigned to the expectation from background. The limited contribution 
of backgrounds to the Z boson candidate sample is reflected in a <1% uncertainty on the final 
measurement. 

The particle-level four-momentum vector of a lepton is computed in the MC simulation by 
adding the four-momentum vectors of any photons found within a radius of AR = 0.1 of each 
lepton axis to the four-momentum vector of the lepton. For the observables of interest in this 
analysis, the use of this cone size makes the electron and muon channels essentially the same 
at the particle level. In this way, the difference in final-state radiation in the Z — ?► mi and 
Z — ?► ee channels is accounted for and the two channels can then be directly combined. The 
particle-level jets in MC events are reconstructed by clustering the generated stable particles 
(after hadronization) using the same anti-Zcj algorithm, with a parameter R = 0.5, as done in 
data. The selection criteria used in data are also applied to particle-level leptons and jets: the 
two leading leptons are required to have px > 20GeV and \ri\ < 2.4, while the jets must have 
Pj > SOGeVand < 2.5. An angular separation of AR (£,;') > 0.4 is also required between the 
two leading leptons and any accepted jet. Finally, the unfolded distributions from the Z — > ee 
and Z ^ jdji channels are combined at the level of covariance matrices using the best linear 
unbiased estimator ||39| . 

The background-subtracted, detector-level distributions are mapped to the particle level by 
correcting for effects of detector resolution and efficiency. Migration of events among bins of 
inclusive jet multiplicity can be caused by detector resolution, especially from the mismeasure- 
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ment of jet pj. For example, an event containing a Z boson produced in association with N jets 
at the particle level, can migrate to the N + 1 jets final state as a result of detector resolution. 
The opposite effect can also occur leading to loss of events that migrate out of the geometric and 
kinematic acceptance. Such migrations correspond to as much as 30% and are treated in the 
detector unfolding procedure summarized below. Detector effects are expressed through a re- 
sponse matrix, which is determined from MC simulation, separately for each lepton flavor and 
each observable, by associating the particle-level values to their reconstructed quantities. Two 
alternative representations of the response matrix, one generated from MadGraph (baseline) 
and the other from SHERPA, are used in this procedure, and half of the difference of the propa- 
gated results is used to define their systematic uncertainty. The unfolding of data to the particle 
level is performed using the Singular Value Decomposition method [40J, implemented in the 
ROOUnfold package [41 J. The total systematic uncertainty due to the unfolding procedure is 
<5% for azimuthal correlations and <2% for the thrust analysis. 

Among the azimuthal observables, the distribution of A(p{Z,ii) has the largest systematic un- 
certainty. This variable is particularly sensitive to the jet-energy scale, which can affect jet 
multiplicity, and A(p{Z,ji) is thereby subject to a larger migration of events from the impact of 
resolution in jet fj. The uncertainty on jet-energy scale varies between 1 and 3% [36], and re- 
sults in an uncertainty of 2-4% on the distribution in A(p{Z, ji), increasing monotonically for de- 
creasing angles. The impact from the resolution on pj is estimated by changing jet resolutions 
by ±10% (corresponding to about one standard deviation) [36J, and comparing the unfolding 
correction before and after these changes. This yields a dependence of ?»1% in the normalized 
azimuthal distributions. The uncertainty from pileup is estimated by changing the cross sec- 
tion for minimum-bias events by ±5%. The resulting uncertainty is 4% for A(p{Z,ii) 0, which 
decreases to a negligible uncertainty for A(p{Z,ii) ^ n. The overall systematic uncertainty on 
A(p{Z,ji) is about 5-6% at values of A(p{Z,ji) ^ and about 2% at values close to n. 

The dominant systematic uncertainty on the thrust distribution, which corresponds to about 
2%, is from the uncertainty in jet-energy scale, and can be understood as follows. When the 
energy scale is increased, more jets enter the two sums in Eq. Q, and both sums tend to shift 
to larger values. Conversely, when the jet-energy scale is decreased, their values decrease. The 
contribution from uncertainty in jet-energy resolution is found to affect the transverse thrust 
by 1%, and the uncertainties from selection efficiencies are <2% for the entire range of Inrj, 
while the uncertainties from pileup and background subtraction have negligible impact. The 
first conclusion is implied in Eq. ([l}, as soft additional pileup energy added to the hard jets 
contributes simultaneously to both the numerator and denominator and, to first order, cancels 
in the ratio. The second conclusion follows from the fact that the transverse thrust is measured 
in the inclusive Z + > 1-jet sample, where the signal purity is almost 99%, and background 
subtraction has therefore only a minimal impact on the measurement of In Tj. 

For pj > 150 GeV the uncertainties on the azimuthal variables and on In tj are evaluated fol- 
lowing the same procedure as described above. However, in addition to the uncertainties orig- 
inating from previously discussed effects, the statistical limitations of the MC samples become 
important and the systematic uncertainty on the result therefore increases. The impact of the 
electron and muon energy scale uncertainties has been assessed and found to be negligible. 

7 Results 

The corrected differential cross sections (normalized to unity) are compared to the predictions 
of MadGraph, SHERPA, powheg Z + 1-jet (NLO), and pythia generators. The differential 
cross sections in the Ai^ and In tj variables are divided by the total Z + jets cross section for the 
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range defined by the lepton and jet kinematic selection criteria, i.e., pj > 20 GeV and \r]\ < 2.4 
for leptons, and px > 50 GeV and \rj\ < 2.5 for jets, and AR > 0.4 for jet-lepton separation. 
The distributions in data and MC are therefore normalized to unity. Figure [s] shows A^{Z,ji) 
as a function of jet multiplicity for inclusive Z + Njets production, with Njets > 1/ 2, and 3. 
Figures |4] and |5] show the A^(Z,;,) and A(p{jj,ji^) distributions, where i,k represent jet indices 
in order of decreasing pj, for Njets > 3- For the sake of comparison, the A(p{Z,ji) results for 
Njets ^ 3 from Fig. |3] are also included in Fig. |4] These distributions characterize essentially 
all the azimuthal correlations in the Njets > 3 inclusive jet-multiplicity bin. Finally, the Z + jets 
distributions in In tj are presented in Fig.|6j For both azimuthal and In tt distributions, the sys- 
tematic uncertainties are presented as bands. The statistical uncertainty from the MadGraph 
MC is displayed as a cross-hatched band for each distribution. 

Overall, the measured distributions in A(p{Z, ji ) agree within uncertainties with the predictions 
from MadGraph. The predictions from SHERPA underestimate the measured distributions by 
about 10% whereas powheg predictions overestimate by about 10%. The disagreements with 
SHERPA and POWHEG (as well as between the two models) become less pronounced at larger 
inclusive jet multiplicities (Fig. |3jl. For Njets = 1/ the Z boson and the accompanying parton 
are completely correlated, and A^(Z,;i) « tt (Fig. [l^). When A(p{Z,ii) <^ n, the presence 
of additional hard QCD radiation is implied. Certain configurations of jets with A(p{Z,ji) < 
tt/2 probe events where the Z boson is in the same hemisphere as the leading jet, and the 
pT of the Z boson is therefore balanced by at least two (or more) subleading jets emitted in 
the opposite hemisphere (Fig. [1J3). The importance of the multiparton LO+PS approach, as 
reflected in MadGraph and SHERPA, can be seen when the data are compared to stand-alone 
PYTHIA at A(p{Z,ii) < 2.5 and Njets > 1- For higher jet multiplicities of Njets > 2 and > 3, 
the distribution in A(p{Z,ii) becomes more isotropic, although a strong correlation remains at 

A(p = TT. 

Within uncertainties, good agreement is observed between the data and MadGraph, SHERPA, 
and POWHEG event generators for Njets > 3. Stand-alone PYTHIA is also consistent with the 
distributions in A(p{Z,i'^) and A(p{i2,i3). In PYTHIA, these high-multiplicity configurations are 
generated exclusively from the PS contribution. The important role of the PS approximation in 
modifying the kinematics predicted from fixed-order calculations is emphasized in POWHEG, 
where its predictive power in a multijet environment (Njets > 3) is evident in Figs. [3|-|5j While 
POWHEG represents an NLO prediction only for the leading jet, and additional radiation is 
modeled exclusively using parton showers, good agreement is observed for data with Njets > 3. 

For the region pj > 150 GeV, the A(p{Z,ji) distributions become more isotropic as jet multi- 
plicity increases. In addition, and contrary to the result for all pj, the angular distributions 
between the subleading jets A(p{ii,j]^) also become isotropic (Fig. [s]?). The improved perfor- 
mance of PYTHIA in this region is consistent with the increased phase space available for par- 
ton emission. A similar observation can be made for distributions in In Tj, which are discussed 
below. The level of agreement between PYTHIA and data for distributions in In tj improves for 
pf > 150 GeV (Fig. |6). 

The corrected normalized differential distributions in In tj for all values of p|, and for pj > 
150 GeV, are displayed in Figs.|6^ and|6]3, respectively, while Figs.|6|: and|6|l provide the corre- 
sponding ratios of data and other models to predictions from MadGraph. The distributions 
for large pj indicate an accumulation of events at values of In tj ~ —2, as could be expected, 
because this region of phase space corresponds to contributions from events with a large spher- 
ical component, corresponding to production of two or more jets. Among the four examined 
models, POWHEG and MadGraph are more consistent with the data, being within 10% of the 
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Figure 3: Normalized A(p{Z,i'[) distributions for the leading jet in the inclusive jet-multiplicity 
bins Njets > 1/ > 2, and > 3: (a) all pf and (b) pf > 150 GeV. Plots in (c) and (d) show 
the corresponding ratios of the data (solid points), and of other MC predictions, relative to 
MadGraph. The ratio for pythia MC is not included in these plots. The error bars on the data 
points represent their statistical uncertainties, the solid yellow shaded band around the points 
represent the sum of statistical and systematic uncertainties taken in quadrature, while the 
cross-hatched (cyan) bands reflect the statistical uncertainties on the MadGraph calculations. 
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Figure 4: Normalized A(p{Z,ji) distributions for the inclusive Njets > 3 jet-multiplicity bin: 
(a) all pj and (b) pj > 150 GeV. Plots in (c) and (d) show the ratios of the data and other MC 
predictions, relative to MadGraph, as described in Fig.lSl 
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Figure 5: Normalized A(p{ji,jj) distributions for the inclusive Njets > 3 jet-multiplicity bin: 
(a) all pj and (b) pj > 150 GeV. Plots in (c) and (d) show the ratios of the data and other MC 
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measured distributions, except at large negative values of In tj for pj > 150 GeV, where ~ 20% 
deviations are observed. The level of agreement of SHERPA with data corresponds to 10-15% 
for most of the bins, while PYTHIA shows discrepancies of > 20%. PYTHIA and SHERPA also 
predict too small values of In tj, especially at values dominated by configurations in which 
the leading jet is produced back-to-back with the Z boson. This yields a larger proportion of 
back-to-back Z + 1-jet events relative to data at small In tj, an effect that can also be observed 
in the A(p{Z,ji) distribution of Fig-lsl 



CMS, ^ = 7TeV,L = 5.0fb 



CMS, ^ = 7TeV,L = 5.0fb 



0.18 
0.16 
0.14 
0.12 
0.1 
0.08 
0.06 
0.04 
0.02 




J 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

- Z/y'^rr, p^>0GeV,N al 

' ^ T jets 


1 1 1 1 _ 


• Data 




- —MadGraph 

ii~i~rf'.r.. 

--SHERPA 1 




' --POWHEG (Z+lj) '"■ 
r PYTHIA 6(Z2) i--!"*" 


I- - — i ~- 


yvi 

- 1 


"i ; 


- . H (a) 
1 




P*l 

tv' 




■t' , 1 , , , 1 , , , 1 , , , 1 , , , 


, , 1 , - 



-12 



-10 



■a 

T3 



0.16 
0.14 
0.12 
0.1 
0.08 
0.06 
0.04 
0.02 




' I ' ' ' I ' ' ' I 
Z/y'^rr, p^>150GeV,N 

• Data 

■■■ MadGraph 

- - SHERPA 
--POWHEG (Z+lj) 
PYTHIA 6 (Z2) 

1- - ■i'^' 



-1 



I 



(b) 



I 



In %^ 



-12 



-10 



In 1^ 



1.4 
1.3 

1.2 
1.1 
1.0 
0.9 
0.8 
0.7 



CMS, ^ = 7TeV, L = 5.0fb ' 
^///, MadGraph stat. uncertainty 



it ^ (c) 



-12 



-10 




I ' ' I 

CMS, fi^ = 7TeV, L = 5.0fb ' 

MadGraph stat. uncertainty 



(d) 



In XT- 



Figure 6: Normalized distributions in InTj for (a) all the Njets > 1 data, and (b) for pj > 
150 GeV and Njets > 1- Plots in (c) and (d) show the ratios of the data and other MC predictions, 
relative to MadGraph, as described in Fig.lSj 



8 Summary 

This Letter reports studies of angular correlations among the objects in Z -i- jets final states. 
The measurements are based on data corresponding to an integrated luminosity of 5.0 fb^^, 
collected with the CMS detector at the LHC in proton-proton collisions at a/s = 7 TeV. 

Azimuthal correlations among the Zboson and the accompanying jets, A(p{Z,ji) and A(p{ji,jj^), 
are measured as functions of inclusive jet multiplicity (Njets > 1/ > 2, and > 3). In addition, 
the transverse thrust event-shape variable In tt is used to characterize the events. Two regions 
of phase space are probed: (i) all events, independent of pj, and (ii) the more highly boosted 
subset of events with pj > 150 GeV. The systematic uncertainties are smaller than those arising 
from statistical sources, which dominate in the extreme regions of phase space. 

The data are compared with predictions from MadGraph, SHERPA, POWHEG Z + 1-jet (at 
NLO), and stand-alone PYTHIA Z -i- 1-jet (at LO). PYTHIA corresponds to the simplest model. 



12 



References 



and is used to gauge the importance of additional corrections from LO and NLO ME formu- 
lations that are interfaced with programs that evolve parton showers. Stand-alone PYTHIA 
provides an adequate description of event topologies when the phase space available for par- 
ton emission is large, e.g., for the highly boosted selection on pj. The MC models that combine 
multiparton QCD LO ME interfaced to parton shower evolution tend to agree with the data. 
The Z + 1-jet ME calculation (at NLO) provided by powheg shows agreement with data at 
large jet multiplicity in the entire phase space probed in this study, despite the fact that, be- 
yond the leading jet, additional radiation comes exclusively from parton showers. 

The measurements presented in this study provide a detailed description of the topological 
structure of Z + jets production that is complementary to existing measurements of rates and 
associated jet multiplicities. As theoretical understanding evolves, these results can be used as 
additional probes of the validity of QCD predictions, while also providing confidence in the 
current MC models as useful tools for the description of SM processes and their application for 
determining background in searches for new phenomena. 
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